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PREFACE 


This  paper  considers  the  problem  of  allocating  weapons  to 
achieve  targeting  objectives  while  simultaneously  minimizing 
aggregate  damage  to  surrounding  nonmilitary  facilities,  each 
of  which  has  an  upper  limit  to  the  damage  it  is  permitted  to 
incur.  A model  is  formulated  which  assumes  only  that  damage  to 
individual  targets  or  associated  facilities  does  not  decrease 
as  the  number  of  allocated  weapons  increases.  An  implicit 
enumeration  algorithm,  based  on  that  of  Lawler  and  Bell  (see 
Reference  [33),  is  described  that  yields  optimal  integer  solu- 
tions. An  example  is  presented. 

This  paper  differs  from  IDA  paper  P-1106  (Reference  [2]) 
in  that  it  presents  the  full  generality  of  the  collateral  damage 
minimizing  model,  whereas  P-1106  describes  a model  (NDM) 
tailored  to  specific  design  requirements.  In  addition,  the  code 
listed  in  the  Appendix  may  prove  a prototype  for  a modified  NDM 
with  greatly  decreased  run  times. 


A TARGETING  MODEL  THAT  MINIMIZES  COLLATERAL  DAMAGE 


One  of  the  assumptions  behind  the  argument  to  employ 
counterforce  targeting  of  strategic  weapons  (the  targeting  of 
an  enemy's  strategic  capability),  as  opposed  to  countervalue 
targeting  (the  objective  of  which  is  the  destruction  of  popu- 
lation and  economy),  is  that  sufficient  destruction  of  strategic 
targets  can  be  achieved  without  causing  appreciable  damage  to 
the  surrounding  nonstrateglc  facilities.  This  paper  presents 
a model  which  addresses  the  following  two  questions:  Given  a 

collection  of  weapons,  potential  almpolnts,  and  a configuration 
of  strategic  targets — each  being  assigned  a minimum  level  of 
damage;  and  nonstrateglc  facilities — each  having  a maximum 
level  of  permissible  damage, 

(A)  Is  there  an  assignment  of  weapons  to  aimpoints  (an 
allocation)  that  satisfies  the  above  two  sets  of 
constraints? 

(B)  Of  all  allocations  satisfying  the  above  two  sets  of 
constraints,  what  is  the  one  (or  a one)  that  minimizes 
the  (perhaps  weighted)  sum  of  the  damage  to  the  non- 
strateglc facilities? 

I.  MATHEMATICAL  FORMULATION 

The  fundamental  elements  of  the  model  are  M strategic 

targets,  henceforth  called  simply  "targets,"  N nonstrateglc 

facilities,  or  "nontargets,"  I different  weapon  types,  and 

J potential  aimpoints  to  which  any  weapon  can  be  directed.  An 

allocation  z is  the  matrix  {z,  ,|l=l,...,I;  j=l,...,J}  where 

z.  , , an  integer,  is  the  number  of  weapons  of  type  1 allocated 
^ J J 

to  aimpolnt  J . 


For  each  target  m,  we  suppose  a real-valued  response 

function  f„(z)  which  represents  the  damage  to  target  m from 
in  ^ 

allocation  z.  We  require  that  f (z)  be  monotonically  non- 
decreasing  in  each  component  of  z,  which  is  an  implicit  assump- 
tion that,  given  any  allocation,  the  allocation  of  additional 
weapons  does  not  result  in  less  damage  to  any  target.  Each 
target  m is  assigned  a real  number,  c^,  which  is  the  minimum 
damage  requirement  (targeting  objective),  i.e.,  for  an  alloca- 
tion z to  be  feasible,  it  must  satisfy  f„(z)  > c , m=l,...,M. 

Similarly,  for  each  nontarget  n there  is  a response  func- 
tion g„(z),  monotonically  nondecreasing  in  each  component  of  z, 
and  a real  number  d^  denoting  the  maximum  damage  permitted  to 
this  nontarget.  Further,  each  nontarget  n is  assigned  a non- 
negative  weight,  or  value,  X^. 

The  nonnegative  integer  w^^  is  the  number  of  weapons  of 
type  i available  for  allocation. 


We  can  now  combine  questions  (A)  and  (B)  into  the  follow- 
ing problem  P: 

N 


Ize  h(z) 

= 

l X g (z) 
^ n n n 

subject  to 

(1) 

n=l 

> 

m=l, . 

..,M  ; 

(2) 

m ^ 

m 

gn^z) 

< 

d 

n 

n=l, . 

..,N  ; 

(3) 

< 

^i 

1=1, . 

T • 

• • j J-  5 

(iO 

^1,J 

e 

1=1,. 

• •jlj  J““l>***>*^  > 

(5) 

Z>^ 

is  the  set  of  nonnegative  Integers.  If  problem  P is 
Infeasible,  then  the  answer  to  question  (A)  is  clearly  "no," 
otherwise  an  answer  to  question  (B)  is  ensured  because  the  num- 
ber of  allocations  which  satisfy  constraints  (4)  and  (5)  is 
finite. 


II.  AN  ALGORITHM 


Problem  P admits  solution  by  Implicit  enumeration.  The 
following  algorithm  Is  based  upon  the  lexicographic  technique 
of  Lawler  and  Bell  (see  Reference  [3]) — though,  unlike  the 
Lawler-Bell  approach,  this  algorithm  does  not  use  binary  vectors. 
We  first  identify  the  matrix  z with  a vector  z.  This  can  be  done 
In  a number  of  ways,  one  of  which  Is  through  the  following 
relationship : 

\ " ^1,J  ^ (J-l)-I;  1»1,...,I;  J=1,...,J  . (6) 

Note  that  this  can  be  reversed  as  follows: 

Zl,j  = Zj^,  1 = k - <(^}*  I,  j * k=l,...,K=I.J. 

where  <x>  Is  the  largest  integer  less  than  or  equal  to  x. 

With  this  In  mind,  we  will  drop  the  circumflex,  and  In  the  dis- 
cussion that  follows,  all  allocations  will  be  vectors  In 
l.e.,  K-dlmenslonal  vectors  with  nonnegative  integer  components. 
We  require  two  binary  relations  between  vectors  In 

Componentwise  (partial)  Ordering: 

We  write  x ^ y — ^k  . . . ,K 

X > y If  X > y and  x,  > y,  for  at  least  one  k. 

Lexicographic  Ordering: 

We  write  x > y If  x,  . > y,  ^ where  k'  = max  (k|x.  ^ y,-}> 

^ h ~ ^ ^ ^ ^ 

and  x>yifx>yorx=y. 

Let  ^ k=l,...,K,  h=k-<(^^)>.  . 

Thus  is  a set  of  allocations  that  satisfy  constraint  (5)  of 
problem  P,  and  clearly  contains  all  allocations  that  satisfy 
constraint  (^),  and  so  must  contain  all  solutions  to  problem  P 
providing  problem  P Is  feasible.  Since  | totally  orders  , 
we  could  emimerate  all  the  points  of  and  find  the  solution 
to  P In  this  manner.  However,  the  monotonicity  of  the 
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objective  and  constraint  functions  will  permit  us  to  skip  over 
many  infeasible  and/or  nonoptimal  points.  To  see  this,  we 
need  some  notation.  Consider  a vector  zz<^.  We  will  denote 
by  z+1  the  vector  x,  if  it  exists,  satisfying 


X e 


S' 


X > z 

0^  r ^ 

Li 

y > z y > X 


(7) 


At  most  one  such  vector  exists,  but  may  fail  to  exist  because 
of  the  boundedness  of  The  vector  z^l  will  be  that  vector 

X,  if  it  exists,  satisfying 

X e (!^ 

Z > X ( ® ) 

~ L 

£ > ^ ^ iS  > Z • 

This  vector  will  always  exist  except  for  z e 0.  The  vector  z* 
will  be  that  x,  if  it  exists,  satisfying 

X E 


X > z 


(9) 


y > X 
~ L ~ 


((y  ^ (y  * £))  = 

Intuitively,  z*  is  the  first  vector  in  ^ following  z (in  the 
lexicographic  ordering)  which  is  not  (componentwise)  greater 
than  or  equal  to  z.  For  some  z,  z*  may  not  exist;  however,  we 
will  adopt  the  following  convention:  For  any  z for  which  z* 

does  not  exist,  we  will  set 

(z«-l)j^  = w^  for  h=\d-(J^y.  I,  k=l,...,K  , 


thereby  ensuring  that  z*^-l  exists  for  every  zeu!^ . Crucial  to  the 

algorithm  is  the  observation  that  for  any  zz'<S' , any  y that 

satisfies  z < y < z*-l  also  satisfies  y ^ z. 

^ L ^ L ^ ^ ^ 


Figure  1 outlines  the  fundamentals  of  the  algorithm.  A 
brief  inspection  of  the  flow  chart  will  make  clear  that  the 
algorithm  must  terminate  after  a finite  number  of  steps . If 
h = ~ upon  termination,  the  problem  is  infeasible,  otherwise 
an  optimum  Integer  allocation  will  always  be  found.  The  order 
in  which  the  constraints  are  examined  was  chosen  because,  for 
certain  applications,  this  order  was  efficient.  However,  we 
make  no  claim  that  this  is,  in  any  sense,  an  optimal  ordering. 

For  other  applications,  a different  sequence  of  constraint 
evaluations  might  well  prove  to  be  better. 

III.  A CLASS  OF  EXAMPLES 

We  will  now  look  at  a class  of  examples  with-  point  targets 
and  nontargets , where  the  destruction  of  any  target  or  nontarget 
is  a binomial  random  variable  with  probability  of  kill  dependent 
on  the  allocation,  but  with  Independent  weapons  effects.  We 
will  use  Cartesian  coordinates  to  specify  location,  in  particular, 
target  coordinates  are  m=l,...,M;  nontarget  coordinates 

are  n=l,...,N;  and  aimpoint  coordinates  are 

J ■!,..., J.  For  response  functions  we  will  employ  "probability 
of  kill"  which  is  computed  as  follows:  Let  p^  . be  the  probabil- 

ity  that  a single  weapon  of  type  1,  allocated  to  aimpoint  J, 
destroys  target  m,  conditioned  on  the  weapon's  arrival  at  its 
destination.  The  probability  that  a type-1  weapon  arrives  at 
its  destination,  its  "reliability,"  is  given  by  p^.  Because  we 
have  assumed  Independence  of  weapon  effects,  it  is  not  difficult 
to  compute  the  total  probability  that  target  m is  destroyed  by 
allocation  z,  which  is 


f (z)  = 1 - ir  ir  (1-p.p^  .) 
" “ 1=1  j=i  ^ 

.n 


'1,J 


Similarly,  we  denote  by  p.  . the  conditional  probability 

1 3 J 

that  a single  type-1  weapon  allocated  to  aimpoint  J destroys 
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START 


7-77-75- l» 


Figure  1.  AN  IMPLICIT  ENUMERATION  ALGORITHM 
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nontarget  n.  Therefore,  the  probability  that  allocation  z 
destroys  nontarget  n Is 


g«(z)  = 1 - 


'n  ~ 


TI  TT  (1-p  p .)  . 

1=1  J=1  ^ 


Although  the  values  of  the  parameters  {p|?  .}  and  {p^  .} 

^ > J 1 J J 

can  be  entirely  arbitrary,  within  the  obvious  limits 


m 


“ i pi.j  i 1 


m=l,...,M;  1=1,...,!;  j=l,...,J  , 


n=l,...,N;  1=1,...,!;  J=1,...,J  , 


we  will  use,  for  tutorial  purposes,  the  following  formulae, 
which  are  not  unreasonable  approximations  to  certain  types  of 
weapon  damage  curves  and  have  been  proposed  by  other  analysts 
(see,  for  example,  Eckler,  Reference  [1],  or  McNolty,  Reference 


1,J 


p?  , = exp)-aj^^jj,[(Xj^-5j)^+(yj^-Cj)^][  m=l,...,M;  1=1,...,!;  J=1,...,J 


= exp  I n=l,...,N;  1=1,...,!;  j=l,...,J 


(10) 


where  all  a,  , 3^  are  nonnegative  real  numbers.  The  param- 
i ,m'  1 ,n 

eters  {a.  } and  {3^  } are  measures  of  the  rate  at  which  weapon 

1 ,m  1 ,n 

effects  decrease  with  distance. 

With  these  conventions,  we  can  now  write  explicitly  the 
problem  P'  which  comprises  this  class  of  examples: 

,N,  and  the  values  of 


Given  nonnegative  weights  n=l. 


CmC[0,l] 

d^c[0,l2 

w^eZ^ 

P.c[0,l] 


m=l, 

n=l, 

1“1, 

1=1, 


,M 

,N 
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a.  > 0 
l,m  — 

®l.n  i 0 


u »v 
n*  n 


1=1, ... ,1;  m=l, . . . ,M 
1=1, ... ,1;  n=l, . . . ,N 
m=l, . . . ,M 
n=l, . . . ,N 
J “1 > • • • > J 


minimize 

z 


m 


subject  to 

Sn<£>  = 

^ ■ 1=1 

J 

Z.  . W.  l“lj  . . . ,I 

J=1  ^ 

Zj  . G 1=~1 , ...  ,1,  J"l,  ...  ,J  . 


m=l,. . .,M 


n=l, . . . ,N 


IV.  COMPUTER  APPLICATIONS 

A FORTRAN  routine  to  solve  problems  of  the  type  given  by 
P'  was  written  for  the  CDC  6400  computer,  and  was  used  to  solve 
the  numerical  example  of  this  section.  (A  listing  of  this  pro- 
gram, along  with  Input  formats  are  given  In  the  Appendix.)  The 
values  of  the  parameters  are  listed  In  Tables  1-6.  The  configura- 
tion of  the  targets,  nontargets  and  almpolnts  is  depicted  In 
Figure  2. 

The  routine  ran  for  five  seconds  to  compute  the  optimal 
solution,  z,  given  In  Table  7. 
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Table  1.  TARGET  PARAMETERS 

M«2 


X 

y 

m 

*^Tn 

m 

1 

0 

.8 

2 

1 

0 

.8 

Table  6.  COMPONENTS  OF  6 


n 


_J 

u 

2 

3 

4 

iTos 

.1 

.1 

o 

VO 

JJ 

GO 

.8 

GO 

CO 

Table  2.  NONTARGET  PARAMETERS 


'’n 

1 

-2 

0 

2 

.3 

2 

-1 

-1 

4 

.3 

n 

3 

1 

1 

6 

• 3 

2 

0 

8 

.3 

Table  3. 

AIMPOINT 

PARAMETERS 

1 

-1 

1 

2 

-1 

0 

J 

3 

0 

0 

1) 

1 

0 

5 

1 

-1 

Table  4.  WEAPON  PARAMETERS 


Table  7.  OPTIMAL  ALLOCATION  £ 

/N 

Z 

J 


1 

CM 

3 

4 

5 

1 

0 

0 

0 

0 

0 

1 

2 

2 

0 

1 

0 

2 

h(z)  « ? 

i.2 

= .28 

f'lU)  ” 

.83 

S^(z)  - .24 

.83 

®3^ 

:£)  - .24 

g^jfz)  -=  .28 

Table  8.  OPTIMAL  ALLOCATION  z' 


’"i 

Pi  

1 

2 

3 

4 

5 

1 

6 

.9 

1 

0 

0 

0 

0 

0 

1 

2 

6 

.7 

2 

0 

0 

3 

0 

0 

Table  5.  i 

COMPON 

1 

1 

lENTS  OF  a h(z')  - 4.5  Ei(z')  = -08 

f,(z')  - .81  g,(z')  » .37 

2 ^ 

1 

.1 

.1 

f2(z')  = .82  g3(z')  = .37 

gj,(z')  = .08 

1 

2 

.5 

.5 
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• TARGET 
O NONTARGET 
□ AIMPOINT 


Figure  2.  CONFIGURATION  OF  THE  EXAMPLE 


It  Is  interesting  to  note  that  if  all  the  are  changed 
to  1.0,  which  is  equivalent  to  removing  the  individual  non- 
target damage  constraints,  then  the  optimal  allocation  is  z', 
given  in  Table  8.  In  this  latter  case,  we  have  reduced  total 
collateral  damage  over  that  given  in  Table  7,  but  only  at  the 
expense  of  considerably  greater  damage  to  two  of  the  nontargets. 
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APPENDIX 


FORTRAN  LISTING  AND  INPUT  SPECIFICATIONS 


FORTRAN  LISTING 


PROGRAM  M0LTM0<INPUT*0UTPUT) 
COHMON/LlHlTS/NUTARtNONON»NOAlMtNOWEAP 
COMM0N/TAR6ETS/XTaR(10) tTTAR(lO) tOEST(lO) 

• COMMON/NOntaR/XMOnuO)  »yN0N(10)  »FACTOR(Io)  fUPNONO(lO) 

COMMON/AiMPNTS/XAlMaoo)  fVAlMJlOO) 
COMMON/E»<EP/l*<NyM(lo)fRELBL(lo>  »EFFTaR(Io»10)  »EFFNON(lot lO) 
COMMON/SCRaTCH/PKT  < 10» 10* 100) tPKN (10*10* 100) * IPRESAL (1000) * 
ITARSURV(lo) *IBTNAL(1000) *BTNts(10) *NTV(lo) *SV* IFLAG.NFLAO*0TNNTS 
2*aTNNTV(lo>tNONFLA6 
CALL  REAOlT 

A Call  calCprb 

^ Call  lexo 

Call  out 
End 
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subroutine  REaoIT 

COHMON/LImITS/NOTaR|NONON,NOAIH*NONE«P 
COMMON/TarOETS/XTaR(10) #YTAR(10) tOEST(lO) 

COMMON/NOnTAR/XNON(10) #YNON(Io) *FACTOR(1o»  tUPNONO(lO) 
COMHON^AlMPNTS^XAlHdOO)  lYAlMdOO) 

COMMON/E«EP/l»»NUMdo>  tRELBLdO)»EFrTAR(10*10)  iEFFnON(10#10) 
COHMON/SCRATCH/PKTdOtlO«lOO),PKNdOtlOf  100)  tlPRESALdOOO)  « 
ITaRSURV(Io)  fIBTNALdOOO)  tBlNTSdO)  fNTV(lo»  fSVtlFLAG,Np|.A6t8TNNTS 
2»0TNNTV  do)  .NONFLAG 
REAOI »NOTaR.NONON,NOA1M.NO«EaP 

I FORMAT(AllO) 

00  5 I»1»n0TAR 

READ10*X1aR(I)  .YTARd)  •OEST(I) 

10  FOrNat(3Fio,6) 

5 Continue 

00  lSIal«NONON 

REA020.XNoN(1) .YNON(I) tFAcTOR(l) .UPNOND(I) 

20  FORNAT<AFiO,6) 

15  CONTINUE 

DO  251-ltNOAlH 
REA030»XAIM(I) .YaIHII) 

30  FORNAT(2flO,6) 

25  CONTINUE 

00100l"i»NOWEAP 
REA040tI"NUHd)  .RELBLd) 

40  FORHAT(I1o*F10.6) 

REA0200.  <EFFTARd,M)  ,M"1,nOTAR) 

REAU300* (EFFNON(I.N) .Nai.NONON) 

200  FORHAT(8F10.6/*2F10.6) 

300  ForMAt(8Fio,6/.2FIO,6) 

100  CONTINUE 
RETURN 
END 
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10 


20 


subroutine  CALCPRB 

COMMON/LlHlTS/NOTAR*NONON«NOAIMtNONCAP 
COHHON/TAROeTS/XTaR(10)  »YTAR(10)  tOESTdO) 

COMMON/NOnTaR/XNON(10»  »YNON(Io) »FACTORtl0) »UPNOND(lO) 
COHMON/AlHPNTS/XAlHdOO)  ,YAIH(100) 

COMMON/E*ER/IWNUM(1o) ,ReUbL(io> ♦efftaR(io»io) »EFFnON(io»1o) 

COHMON/SCRaTCH/PKT(10»10»100) tPKN(lOtlO«lOO) tiPRESALdOOO) * 
1TaRSURV(Io{ »IBTNAL(1000) *BTNts(10) *NTV<1o) *SVtXFUA0»NFLA6.BTNNTS 
2*BTNNTVUo*  «nonflao 
0010Mm1«NoTaR 
ooioi«i*nomeap 

OOlOJalt'^OAlM 

NK«EFFTAR(I»M>«( (xAIM(J)-xTAR(N) )*»2*(YAiM<J)-YTAr(m) J««2) 

PKT(NiItd)»R£LBUl)»EXP(-WlO 

CONTINUE 

OOBONaltNoNON 

OO20laltNOWEAP 

Oo20jal«NoAIH 

WW"EFFN0N(I,N)«( (XAIM(J)-XNON(N))««2*(YAiN(J)-YNON<N) )««2) 

PKN (N» I »d» ■RELBL ( I ) »EXP (-WN) 

continue 

Return 

End 
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subroutine  lEXO 

This  subroutine  finds  anu  stores  the  optimal  allocation  using 
lexicographic  enumeration  after  LAWLER*tELL.  OPTIMAL  VALUES  ARE 
STORED  AS  FOLLOWS— 

iTNNTS-«TOTAL  NONTAROET  SURVIVAL  LEVEL  (FROM  NTS) 

■999999999.  IF  NO  FEASIBLE  SOLUTION  Is  FOUND 
iBTNALt.*.)— OPTIMAL  ALLOCATION 
bTnts(.I>«reSUltino  Target  survival  level  (from  TaRS) 

8TNNTV(.)— individual  NONTaRGeT  survival  levels  (FROM  NTS) 
COMMON/LlMiTS/NOTAR.NONON.NOAlMtNOWEAP 
COMMON/tAROET$/XTAR(IO)fYTAR(10) .OESTIlO) 

COMMON/NOnTAR/XnOnUO)  tVNON(lo)  fFACTOR(lo)  tUPNONO(lO> 
COMMON/AiMPNTS/XAlM(lOO),TAIM(lOO) 

COMMON/E»EP/IWNUM(lo).RELBL(lo)fEFFTAR(lo»IO)»EFFNON(ioflO) 
COMMON/SLraTCH/PKT ( 1 0 1 10 1 1 00 ) «PKN ( 10 1 1 0 ( 1 00 ) t IPRESAL  (1000). 
lTARSURV(lol.lBTNAL(1000)tBTNTS(10) .NTV(lo>  *SVtIFLA6.NE|.AG*BTNNTS 
2«BTNNTV(1o) tNONFLAO 
REAL  NTV 

INTEGER  ITCMPAL(IOOO) 

INTEGER  ITEMST(IOOO) 

00  9000  Ll«1.M 
ITEMST(LL)*o 
9000  CONTINUE 

BTNNTS»9999999999. 

MbNOVEaR*NOaIM 
DO  1 K*ltM 
IPRESAL(K)bo 
1 CONTINUE 

C BEGIN  enumeration 

C CHECK  EERo  vector  FOR  FEASIBILITY 

CALL  TARS 

IF(IFLAG.nE.0)OO  to  100 

C IF  HEREt  NO  ADDITIONAL  VEaPONS  aRE  NEEDEq 

8®0E  call  NTS 

IF (NONFLAG.NE.o)RETURN 
8003  BTNNTSaSV 

DO  10  KMl.NONON 
BTNNTV(K>bNTV(K) 

10  continue 

00  11  K«1,M 
1BTNAL(kU1PRESAL(K) 

11  CONTINUE 

00  12  KmI.NOTAR 
8TNTS(K)«TARSURV(k) 

12  continue 
^cturn 

c This  section  computes  next  allocation 

3l0  DO  315  U"1(N0AIM 
00  315  i"l.NOVEAP 
KKK»1*(J-1)«N0wEAP 

1F(IPRESAl  (KKK).LT.IWNUM(I))oO  to  320 

IPRESaL(Kj<K)«0  

315  CONTINUE 

C here  if  IpRESAL  WAS  LAST  ALLOCATION 

Return 

320  IPR£SaL(KkK)*IPRESAL(KKK).1 
399  CALL  NUMBS 


IF(NFLA6.CQ.1)00  to  600 

*00  Call  tars 

IF(IFLA0.nC,1)O0  to  500 
C here  if  allocation  infeasible 

C STORE  IPftESAL 
loo  Do  *05 

ITEMPAL(^)■IPRESAL(K) 

*05  continue 

C NOW  TO  compute  IPRESALSTAR-1 
00  *10  H«itM 

lF(IPRE§AL(lO  .E0»0)00  To  *l0 
IPRESAL(K).o 
00  TO  *15 
*10  CONTINUE 
*lS  IFIK.GE«H)00  TO  420 
LbK^I 

00  425  F*LtM 

J«(K-1)/N0WEAP 

l«K»NOwEAp»j 

IF(1PRESAl(K).LT.IWnUM(I) )00  TO  430 
1PRESaL'(K)«0 
425  CONTINUE 
GO  TO  420 

430  lPRESAL(K)alPRESAL(K)«l 
00  TO  435 

*20  DO  **0  l"ltNOMEAP 
DO  **0  J"1»N0aIM 
KKK»1^( j-l>«NOKEAP 
IPRESALIKkK) alUNUM ( 1) 

**0  CONTINUE 
00  TO  *&0 

*35  00  **5  j"lfNOAlM 
DO  **5  I"1#N0WeaP 
KKKaU(J-l)«NOMEAP 
1F(1PRESAl(KKKJ *NE.0)OO  TO  *50 
IPRESaL<KkK)»1WNUm(I) 

**5  CONTINUE 

*50  lPRESALlKKK)aIPRESAL(KKK)-l 
C NOW  WE  HAVE  IPRESaLSTaR  -1 

46o  do  9005  LL-ltH 

If (1TeMST(LL>*NE«IPRESaL(LL))gO  to  9010 
9005  continue 
60  TO  902o 
9010  DO  9015  LLaitH 

ITEMST(LL)"1PReSaL(LL) 

9015  continue 

Call  taRs 

IF(IFLA0.nE.0)60  to  310 
C here  IFIpReSALSTaR  -1  IS  feasible 

9020  CONTINUE 

DO  *95  F"1»M 
IPRESaL  <N) "ITEMPAL (k) 

*95  continue 
GO  TO  310 
500  CALL  NTS 

IF (NONFLA6*NE*0)OOT0600 
BOlO  IF(SV,OE«BTNNTS)00  to  600 

c Have  foONq  a ne*  optimum 
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c 


00  510  K»i»M 
1bTNaL<R>»1PRESaL(K) 

510  continue 

00  515  I"1»N0TAR 
8TNTSn)*TARSURV(i) 

515  CONTINUE 
BTNNTSaSV 
00  520  i"l(NONON 
BTNNTVa)«NTV(l) 

520  continue 

Skip  to  Ipresalstar 

600  00  610  K»1(H 

1P(IPRESA|,(k)  .EQ»0>00  to  6^0 
IPRESAL(K)-o 
QO  TO  620 
610  CONTINUE 
620  IF(K.GE«H)  return 
L»K*1 

00  625  KaL»H 

J«(K-1)/NoKEAP 

IbKoNOmEApaj 

IF(IPRESAl(K» .LT.IKnUM(i) )Q0  TO  630 
lPRESAL(K)aO 
625  CONTINUE 
RETURN 

630  IPRESaUK)»IPRESAL(K)  *1 
00  TO  399 
End 


A-6 


subroutine  tars 

COMHON/LlMiTS/NOTARtNONON.NOAIHtNOWEAP 
COMHON/TaROETS/XTaR(10)  tYTARdO)  iDEST<10) 

COMMON/NOnTAR/KNON (10) » YNON ( 1 0 ) t FACTOR ( lo ) »UPNOnO ( 1 0 > 
COMHON/AiMPNTS/XAlHdOO)  tYAIM(lOO) 

C0HM0N/EneP/IMNUM(10)  tREL8LdO)*EFFTAR(lOtl0)  lEFFNONdOi  10) 
COMHON/StRATCH/PKTdOt  10(100)  tPKNdOdOtlOO)  dPRESALdOOO) , 
ITARSURV(Io) (iBTNALdOOO) (BTNts(IO) (NTV(Io) f SVi IFLA6(NflA6»BTNNTS 
Z(BTNNTVdO)  (NONFLaG 
OOlOMalfNoTAR 
TSURVal* 

OOSOlal*NOMEAP 
DOSOjaltNoAIH 
KKKai*(J»l)*NO)(eAP 
lFdPRE?AL(KKK)  .LE»^>GOT050 
IF  (PKT(M(ltj)  *0Ed»)B(9 
a TSURVaO* 

GO  TO  SO 

9 PSa(l,-PKT(Md,J))«»IPRESAt(KKK) 

TSURVaTSURVapS 

SO  CONTINUE 

C ARE  constraints  SATISFIED 

WMal,-OEST(H) 

1F(TSURV.gT,RW)6»7 

6 iFLAOal 

return 

7 TARSURV(M)«1,-TSURV 

10  CONTINUE 
IFLAGaO 
RETURN 
End 
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subroutine  NTS 

COMMON/Ll.HlTS/NOTAR*NONON,NOAlMfNOWCAP 
COHMON/TAROETS/XTaR(10) »YTAR(10) tOEST(lO) 

COMMON/N0nTaR/xN0N(10) »YNON(io) »FACTOR(lo) »UPNOND(lO) 
COHHON/AimPNTS/XAIHUOO)  tYAIHaQO) 

COMMON/£PeP/1NNUM(1o)»RELBL(1o)»CFFTAR(1o»1o>  »EFFNON(io»10) 
COMMON/SURaTCH/PKT(10»10»100) •PKNdOtlOtlOO) •IPRESAL(IOOO) * 
ITARSURV(IO) *IBTNAl(1000) ♦PTNts(Io) »NTV(lo) *SV»1FLaG»NflA6»BTNNTS 
2*BTNNTV(1o) »nonflao 
REAL  Ntv 
SVaO* 

OOlOOONal *NONON 
TEHPSVal* 

00501"1»^OWEAP 
OOSOJalfNOAIH 
KKKaI«(J-l)*NOMEAP 
IF<IPRESAl(KKK) .LE«o)GOT050 
lF(PKN(N*i,j) .GE«1.)8*9 

8 TEMPSVaO, 

GO  TO  SO 

9 PSa(l,-PKN(N»I*J) )»«IPRESaL(KKK) 

TEMPSVaTEMPSVaPS 

50  CONTINUE 

WNal.-TEMpSV 
IF(KW.GT»UPN0ND<N) )7»10 
7 NONFLAOal 

Return 

10  SVaSV^FACTOR<N)aWH 
NTV(N) aFAcTOR(N)»RR 
CONTINUE 
NONFLAGaO 

RETURN 

End 
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subroutine  numbs 

COHHON/LlHlTS/NOTAR*NONON«NOAlMtNOMEAP 
COMMON/TaRQe]S/XTaR(10)  ♦YTARUO)  .OEST(IO) 

COMMON/NOnTAR/XNON(10) tYNON(lo)#FACTOR(lo) tUPNONO(lO) 
COMMON/AiMPNTS/XaIM(100),Ya1M(100) 

COMMON/E«eP/ I «NUM ( 1 0 ) t ReLbL ( 1 0 ) • EFFTaR ( 1 0 ♦ 1 0 > • EFFnON ( 1 0 » 1 0 > 
COHMON/ScRATCH/PKTdOtlOtlOO) tRKN(l0(l0*100) tlPRESAL(lOOO)* 
iTARSURV(lo)»lBTNAL(IOOO) »BTNts(10) iNTV(lo) »SV. IFLaG*Nfl*G,BTNNTS 
2t8TNNTV(lo)*NONFLAG 
0011»1»N0wEaP 
ISUH>0 

OOZ^iiltNOAlM 

KKK*l4<J-l»»NOWEAP 

isuhmipResal (KKK) ♦ISUH 
continue' 

IF<ISUh;Le*INNUM(i) )60T01 

nflag-1 

RETURN 

CONTINUE 

Nn.AQ«0 

RETURN 

END 


10 

15 

20 


25 

26 

111 

30 

35 


36 

4l 


45 

SO 


55 

56 


60 


65 


70 

71 

100 


101 


subroutine  out 

COHHON/UMlTS/NOTARtNONON»NOAlHfNOMEAP 
C0MM0N/TAR0ETS/XTaR<10) lYTAR(lO) tOEST(lO) 

COMHON/NOnTaH/XNON(10) »YNON(1o) tFACTOR(lo) »UPNONO(lO) 
COHHON/AiHPNTS/XAlMdOO)  tYAlM(lOO) 

COMMON/EPEP/IWNUM(1o) fRELBL(lo) tEFFTAR(lo»lO) »EFFNON(lotlO) 
COHHON/SCRATCH/PKT(10»10ilOO)iPKN(10tlOtlOO) »IPrEsAL (1000) « 
iTARSURVUAltlSTNALdOOOltBTNTSdo)  tNTV(lo>*SVtIFLA6«NFLA6.8TNNTS 
2»bTNNTV(*0) .nonflao 
Real  ntv 
PR1NT5 

format (IHi.«PROORaM  mol*) 

PRINTlO 

FORMAT(lM-t*INPUT  DATA*) 

PRiNTlStNoTAP 

F0RMAT(1«-#«TaR0ETS  (•tI4,*)«) 

PRINT20 

format (lM()f«TARGET  NUM8EH«»6X»*X  COORO.«t7X,«Y  COORO»«f6Xt 
^•PROB.  OF  0EST.«/) 

0026I»1»MaR 

PRINT25f IfXTAR(I) .YTaR < I ) .OEST ( I ) 

FORMAT (IH  t7Xfl2tl0XfFl0.3f5X,Fl0.3t5X«F10*3) 

CONTINUE 
PRINTlll 
format (iM-) 

PRINT3O1N0NON 

FORMAT  dH«#«N0NTAR6ETS  (•♦14#*) •) 

PR1NT35 

Format (lMof*NONTAR  NUM8EH*t6X#«X  COORD. •»7X,*Y  COORO.«»12X. 
mvaLUe**  7X. •Damage  limit*/) 

D0411>lfN0N0N 

PRINT36*I«XN0N(I) yYNON(I) .FACTOR(I) .UPNOnDCD 

FORMAT  (iM  »7X,I2.ioX#Flo»3»5X.Flo«3.5X,FiOOf5X.Fio.3) 

CONTINUE 
PRINTlll 
PRINT45. NOAIM 

FORMAT(lM.t«AlMPOlNTS  (•tI4*)#) 

PRINT50 

FORMAT (lHOt«AlMPNT  NUMaEH*t6X,*X  COORD. •y7Xy*Y  COoRD**/) 

D056l»l »N0AIM 

PRINT55#I.XaIM(I) .Y*IM(i) 

format (IH  t6XtI3f 10XtFl0*3t5X.Fl0.3) 

Continue 

PRINTlll 

PRINT60»NOMEAP 

FORMATdn-. •WEAPON  CLASSES  (•.I4.^)^) 

PRINT65 

FORMAT (lMo,*CLASS  NUMBER^,5X,*TOTAL  AVAILABLE*. 5X, •RELIABILITY*, 
l/> 

0071I*1«N0WEAP 
PRiNTTO.I.IWNUHd)  tRELsL  (D 
FORMAT (IM  .5X,l2,l5X,I4,l5X,F10.3) 
continue 
PRINTIOO 

format (IM-»*WEAPOn-TAROET  EFFECTIVENESS  TABLED) 
print  lol 

FORMATdM-.*  /TARGET*) 
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PR1NT102 

102  FORHATilH  t«V<EAPON/«) 

PRINTl03t(I»I*l»NOTAR) 

103  rORMAT(lHOfl^XfI2,9<8X.I2)/) 

OOl051"ltNOWEAP 

PRINTUQ»lt(EFFTAR(ItJ)»J*lfNOTAR) 

UO  F0RHAT(1«  tAXfl2f4X,l0Fl0,3) 

105  continue" 

PRINTlli 

PRINT200 

200  FORMAT (lM-,»t(EAPON-NONTAHGET  EFFECTIVENESS  TABLE*) 

PR1NT201 

201  FORMAT (lM-,«  /nONTAHGET*) 

PRINT  102 

„ PR1NT2o3» (I.I«1»N0N0N) 

203  FORMAT<lM0|l4X»l2t9(8X»l2)/) 

Oo2051m1«n0wEAP 

PRINT210»I,(EFFNON(1.J) *J«1,N0N0N) 

2l0  F0RMAt<1M  t4X*I2tAX*10Fl0.3) 

2U5  CONTINUE 

PRINTUI 

PRlNTin 

PRINTIOOS 

1005  format (iH-f •allocation  RESULTS*) 
lF(8TNNTS,tT.9999999999.)60TOllOO 
PRINTlOlO 

1010  formaT(iMo,«it  is  Impossible  to  meet  the  target  damage  constraints 
1— problem  infeasible*) 

Return 

1100  PRINT  1125 

1125  FORMAT (lMoi*FOLLOWlN6  IS  TME  OPTIMAL  ALLOCATION#) 

PRINT1130 

1130  Format (lM0f«WEAPON  CLASS*»5X»«AlMP0lNT»t5X*«NUMBER  assigned*/) 
D01l36I«ltN0WEAP 
001136J*1*N0AIM 
KKK»I» (J»1)*N0WEAP 
IF<I8TNaL<KkR).LE.O)GOTOH36 
PRlNT1135tIt J(I8Tnal<KKK) 

1135  FOHMATdH  *5X  • l2«  l3x*  l3,  loX  1 1 10) 

1136  continue 
PRINTlli 
PR1NT1205 

1205  FORMAT (1H.*«TARGET  DAMAGE*) 

PRINT1210 

I2l0  F0RMAT<1N-»*TAR0ET  CLASS*»5Xt«REsULTlNG  PROB.  OF  0EsT,«i5Xt*SPEClF 
llEO  PROB*  of  OEST,*/) 

00  l2l6l*itN0TAR 

PRINT  I215»|,(jTnTS<I)  *oe5T(D 

1215  foRmATUH  tZXtj2»i4XtFl0.3»17X*Fl0.3> 

1216  continue 

PRINTlli 

TOT"0» 

OO12201*liNONON 

TOT»TOT*FaCTOR(I) 

1220  CONTINUE 

PRINTl225,T0T 

1225  FORMAT<lHO»*ORIOlNAL  TOTAL  NOnTARGET  VaLUE  WAS  *»f10;3) 
PC"BTNNTS/TOT*100, 


PR1NTi230*BTNNTS»PC 

1230  format (1M-,*T0TAL  EXPEcTeo  NONTaRGET  Value  destroyed  is  ♦.Flo. 3. 

1*  or  •|F16.3.*  percent.*) 

PRINTUI 

PRINT1235 

1235  Format (iH-*«iNoiviooAL  nonTarget  expected  Value  oesthoyed  listed  8 
lELOV*) 

PRINT1240 

1240  F0RMAT(IH-,»N0NTARGET  NuMbEr*,5X.*0R1GINaL  value*. 5Xt 

l*EXPECTEO  Value  destroyed*. 10X.*PERCENT*,5X.*SPEC1FIEd  maximum  (PE 
2RCENT)*/) 

DO  I246lal.N0N0N 
4WW»1oO.*UPNOND(I) 

PCC*BTNNTv(I)/FACTOR(I)*1oO. 


1246  CONTINUE 


Return 

END 


INPUT  SPECIFICATIONS 


Refer  to  problem  P"  for  notation. 

Core  requirements  Impose  the  following  limits: 


Input  Parameters  Format 


M < 10 
N < 10 
J < 100 
I < 10  . 

Card  Name 

LIMITS 

TARGET 

(1  card  each) 

NONTARGET 
(1  card  each) 

AIMPOINT 
(1  card  each) 

WEAPON 

(1  deck  each) 


M,  N,  J,  I 


^m> 

^ J 

X , d 
n n 

'j 

“l,l’  “l,2’ 
^i,l’  ^i,2’ 


illlO 

3F10.6 

4F10.6 

2F10.6 

I10,F10.6 

8F10.6/2F10.6 

8F10.6/2F10.6 


• • • 


